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Abstract — Recent results establish the optimality of interfer- 
ence alignment to approach the Shannon capacity of inter- 
ference networks at high SNR. However, the extent to which 
interference can be aligned over a finite number of signalling 
dimensions remains unknown. Another important concern for 
interference alignment schemes is the requirement of global 
channel knowledge. In this work we provide examples of iterative 
algorithms that utilize the reciprocity of wireless networks to 
achieve interference alignment with only local channel knowledge 
at each node. These algorithms also provide numerical insights 
into the feasibility of interference alignment that are not yet 
available in theory. 



I. Introduction 

The recent emergence of the idea of interference alignment 
for wireless networks has shown that the capacity of wireless 
networks can be much higher than previously believed [ 1 ] . The 
canonical example of interference alignment is a communica- 
tion scenario where, regardless of the number of interferers, 
every user is able to access one half of the spectrum free from 
interference from other users [ 1 ] . For the interference channel 
with K transmitters and K receivers and random, time varying 
channel coefficients drawn from a continuous distribution, 
reference [1] characterizes the network sum capacity as 

C S {SNR) = — log(SNR) + o(log(SNR)) (1) 

so that the capacity per user is | log(SNR) + o(log(SNR)). 
Here SNR is defined as the total transmit power of all the 
transmitters in the network when the local noise power at 
each node is normalized to unity. The o(\og(SNR)) term, by 
definition, becomes negligible compared to \og(SNR) at high 
SNR. Therefore the accuracy of the capacity approximation 
in (1) approaches 100% at high SNR. Since the capacity of a 
single user in the absence of all interference is log(SNR) + 
o(\og(SNR)), the main result of [1] may be summarized as: 

"Af high SNR, every user in a wireless interference network 
is (simultaneously and almost surely) able to achieve (nearly) 
one half of the capacity that he could achieve in the absence 
of all interference." 

The capacity-optimal achievable scheme within 
o(log(SNR)) is shown to be the interference alignment 
scheme. Interference alignment on the K user interference 
channel refers to the idea of constructing signals in such a 
way that they cast overlapping shadows over one half of the 



signal space observed by each receiver where they constitute 
interference, leaving the other half of the signal space free of 
interference for the desired signal. This approach reveals the 
sub-optimality of the cake-cutting view of spectrum allocation 
between co-existing wireless systems because, essentially, 
everyone gets "half the cake". 

Interference alignment schemes are presented in [1] in the 
form of closed form expressions for the transmit precoding 
matrices. However, these closed form expressions require 
global channel knowledge which can be an overwhelming 
overhead in practice. Moreover, closed form solutions have 
only been found in certain cases. In general, analytical solu- 
tions to interference alignment problem are difficult to obtain 
and even the feasibility of interference alignment over a 
limited number of signalling dimensions is an open problem. 
In this paper we explore distributed interference alignment 
algorithms to accomplish the following objectives. 

• Require only local channel knowledge at each node. 
Specifically, each receiver is assumed to know only the 
channel to its desired transmitter and the covariance 
matrix of its effective noise (consisting of the AWGN 
and the interference from all other users). 

• Provide numerical insights into the feasibility of align- 
ment. 

We propose iterative algorithms that take a cognitive ap- 
proach to interference management and utilize only the local 
side information available naturally due to the reciprocity of 
wireless networks. The two key properties can be summarized 
as follows. 

• Cognitive Principle: Unlike selfish approaches studied 
in prior work where each transmitter tries to maximize 
his own rate by transmitting along those signaling di- 
mensions where his desired receiver sees the least in- 
terference, we follow an unselfish approach where each 
transmitter primarily tries to minimize the interference 
to unintended receivers. Since avoiding interference to 
unintended receivers is the defining feature of cognitive 
radio [2] the unselfish approach is a cognitive approach. 
The cognitive approach is found to lead to interference 
alignment, and is thus capable of approaching network 
capacity at high SNR. 

• Reciprocity: For a given transmitter, learning how much 
interference is caused at unintended receivers can re- 
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quire too much side information, and is one of the 
key challenges for cognitive radio systems. However, 
this information is naturally available because of the 
reciprocity of the channel for networks where two-way 
communication is based on time-division duplex opera- 
tion with synchronized time-slots. Due to reciprocity, the 
signalling dimensions along which a receiving node sees 
the least interference from other users are also the same 
signalling dimensions along which this node will cause 
the least interference to other nodes in the reciprocal 
network where all transmitters and receivers switch roles. 
The paper is organized as follows. In the next section, 
we review some optimization approaches for interference 
networks in existing literature. In Section IV we state the 
open problem of determining the feasibility of interference 
alignment over a limited number of signaling dimensions. 
The same section also describes the reciprocity property of 
interference alignment. In Section V we present iterative 
interference alignment algorithms. In Section VI, we show 
that these algorithm can achieve performance close to the 
theoretical results and discuss a few applications. We conclude 
with Section VII. 

Notation: We use lower case for scalars, upper case for 
vectors and bold font to denote matrices. represents the 
d th column of matrix A. 1^ represents the d x d identity 
matrix. Tr[A] denotes the trace of the matrix A and A^ is the 
conjugate transpose of matrix A. Finally, K, = {1, 2, • • • , K} 
is the index set of K users. 

II. Interference Optimization Approaches 

The optimality of interference alignment schemes at high 
SNR is interesting because these schemes treat all interference 
as noise and require no multi-user detection. Achievable 
schemes based on treating interference as noise have been 
explored extensively over the last decade. Prominent among 
these are the interference avoidance and iterative waterfilling 
algorithms where each transmitter acts selfishly to align its 
transmissions along those directions where its desired receiver 
sees the least interference [3]-[7], and network duality ap- 
proaches [8]— [ 11] that are based on the reciprocity of the 
wireless propagation channel. 

A. Interference Avoidance and Iterative Waterfilling 

Iterative algorithms are commonly used for various resource 
allocation problems, such as interference avoidance and it- 
erative waterfilling. However, the philosophy of interference 
alignment is quite distinct from both iterative waterfilling and 
interference avoidance. With iterative waterfilling/interference 
avoidance algorithms [3], [12], [13], each transmitter tries to 
do what is best for his own receiver, i.e., each transmitter 
allocates its power in a manner best suited for his desired 
receiver. With interference alignment each transmitter tries 
to minimize the interference he causes to other receivers. 
The interference alignment schemes in [1], [14]— [16] show 
that for interference networks, the "do no harm" approach is 
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Fig. 1 . Interference alignment solution for the three user two antenna case. 
The arrows in the red indicate the direction of the interference. 

much more powerful, and is in fact capacity-optimal within 
o(log(SNR)), than the "help yourself" approach of interfer- 
ence avoidance and iterative waterfilling schemes. 

From a game-theoretic perspective, interference avoidance 
and iterative waterfilling algorithms lead to a stable operating 
point commonly known as the Nash equilibrium. At Nash 
equilibrium, there is no incentive for any user to unilaterally 
change his transmit strategy. Often these points are not op- 
timum from a network perspective and indicate inefficiency 
in wireless network operation [17], [18]. Interestingly, inter- 
ference alignment is not a Nash equilibrium point if the goal 
of each user is to maximize his own rate. Fig. 1 shows the 
interference alignment solution for the three user two antenna 
case. Notice that interfering signals are co-linear at each 
receiver while the desired signal may not be exactly orthogonal 
to the interference - a price paid for interference alignment. It 
can be easily observed that the interference alignment solution 
is not a Nash equilibrium point. Fixing the transmit strategy 
for users 2 and 3, the best strategy for user 1 is to choose 
Ui such that his signal is orthogonal to the interference at 
receiver 1. Although this strategy is good for user 1, it will 
destroy interference alignment at receivers 2 and 3. Thus Fig. 
1 clearly highlights the difference between the optimal strategy 
of interference alignment and the selfish strategy of iterative 
waterfilling or interference avoidance. 

Iterative schemes have also been used to implicitly achieve 
interference alignment on the 2 user X channel in [19]— [21]. 
However, for the 2 user X channel interference alignment 
can be explicitly achieved with roughly the same amount 
of channel knowledge as required by the iterative schemes, 
without the need for an iterative process [14]. The iterative 
schemes of [19] are specialized for the 2 user X channel and 
generalizations to X networks and interference networks with 
more than 2 users are not straightforward. 

B. Network Duality 

Another approach taken in prior work is to exploit the dual- 
ity relationships enabled by the reciprocity of the propagation 
channel. For example, network duality ensures that the same 
set of signal to interference and noise ratios (SINRs) can be 
achieved in the original and the reciprocal network with the 
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same total transmit power [8], [9]. Network duality is used 
in [9], [10] to minimize the total transmit power required 
to support a feasible rate vector. Reciprocity of propagation 
channels is used in [11] for optimal frequency allocation 
problem. 

In this work we provide examples of iterative algorithms 
to achieve interference alignment on wireless interference 
channels. These algorithms combine elements of all the above- 
mentioned approaches, especially [19] and [11]. 

III. System Model 

Consider the K-usei MIMO interference channel where the 
k th transmitter and receiver are equipped with and 
antennas respectively. Note that the antennas could represent 
symbol extensions in time or frequency as well. However, if 
the antennas correspond to symbol extensions over orthogonal 
dimensions (time, frequency slots) then the channel matrices 
will have a diagonal structure. The channel is defined as: 

K 

yW(n) =^H' H l(n)l' ! l(n) + Z [fe] (n), Vfc e K 
i=i 

where, at the n th channel use, (n), Z^ (n) are the N^xl 
received signal vector and the zero mean unit variance circu- 
larly symmetric additive white Gaussian noise vector (AWGN) 
at receiver fc, X^(n) is the M' 1 ' x 1 signal vector transmitted 
by transmitter I, and H^ kl '(n) is the iVM x MM matrix of 
channel coefficients between transmitter I and receiver k. The 
transmit power at transmitter / is E[||X W || 2 ] = P [l] . 

For the K user interference channel defined above, we also 
define a reciprocal channel, where the role of transmitters 
and receivers are switched. For every variable on the original 
channel, the corresponding variable on the reciprocal channel 
is denoted with a left arrow on top. The reciprocal channel is 
defined as: 

K 

Y I fe l (n) = (n) + *Z [k] (n), Vfc e K 

i=i 

where, at the n th channel use, Y^{n), Z^(n) are the 
M^xl received signal vector and the zero mean unit variance 
circularly symmetric additive white Gaussian noise vector 
^AWGN) at receiver k which is equipped with antennas, 
X™(n) is the x 1 signal vector transmitted by transmitter 
I, and H[ fe 'l(n) = U^(n) is the x matrix of 

channel coefficients between transmitter I and receiver k. The 
transmit powers_at transmitter I on the reciprocal channel is 
E[||JfM|| 2 = PW. The channel use index n is henceforth 
suppressed to avoid cumbersome notation. 

IV. Interference Alignment over Limited 
Dimensions - An Open Problem 

References [1] and [15] present interference alignment 
schemes that are constructed over symbol extensions of time- 
varying channels. It is shown that by using long symbol exten- 
sions the degrees of freedom achieved per dimension approach 



arbitrarily close to the theoretical outerbound, thereby estab- 
lishing the degrees of freedom of time-varying interference 
and X networks. However, the extent to which interference 
can be aligned over a limited number of dimensions remains 
an open problem. As a consequence, the maximum number 
of degrees of freedom that can be achieved through alignment 
of interference signal vectors is not known in general. Note 
that interference alignment can also be accomplished in terms 
of signal levels rather than signal vectors by using structured 
codes (multilevel and lattice codes) as shown in [22], [23]. 
However, in this work our focus is on interference align- 
ment through signal vectors. We first review the interference 
alignment problem and the reciprocity property of interference 
alignment. 

Let $® < min(M[ fc ], N^),k e JC denote the degrees of 
freedom for user fc's message. 

Precoding at Transmitter: Let V[ fe ] be an x d\ k ^ matrix 
whose columns are the orthonormal basis of the transmitted 
signal space of user fc. Mathematically, the transmitted signal 
vector of user fc is given by: 




Each element of the d} k ^ x 1 vector X represents an 
independently encoded Gaussian codebook symbol with power 
^pq- that is beamformed with the corresponding vector of . 

Remark: For interference alignment or the achievability of 
the degrees of freedom it suffices if the beamforming vectors 
are linearly independent. However, we assume the beamform- 
ing vectors are orthonormal in the formulation above. Note that 
this does not affect the feasibility of interference alignment, 
and it naturally leads to the iterative algorithm to be presented 
in this paper. 

Interference Suppression at Receiver: Let be an x 
c?[ fe ] matrix whose columns are the orthonormal basis of the 
interference-free desired signal subspace at receiver fc. The k th 
receiver filters its received signal to obtain: 

Y [k] = u [fe]t r [fel (3) 

If interference is aligned into the null space of Tjl fc l then the 
following condition must be satisfied: 

uWtH^V^ = O.Vj^fc (4) 
rank^urotH^vW) - d k (5) 

In other words the desired signals are received through a d^ x 
dl k ) f u u rank channel matrix 

jjlfefe] a u[fc]t H [fcfc] v W 



while the interference is completely eliminated. The effective 
channel for user fc is then expressed as: 

where Z W - JV (0, 1) is the effective #1 x 1 AWGN vector 
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at receiver k. The rate achieved on this channel is: 



Rlk] 



log 



P [k] W [fcfc]Tf[fcfc]t 



(7) 
(8) 



= d[ fe llog(p[ fe ])+o(log(p[ fc ])) 
Thus, d-W de grees of freedom are achieved by user k. 

A. Feasibility of Alignment 

Given the channel matrices H^l, k, j G JC, we say that the 
degrees of freedom allocation (dM , <P^ , • • • , ) is feasible 
if there exist transmit precoding matrices and receive 

interference suppression matrices U' fe ': 



such that 



tjM . j^W x d M 

U[fe]tjj[fcj] y[j] 

rank (u^H^ V^) 



V [feIt V [fe] = I dW (9) 

UlfcltuW = IrfW ( io) 

0,Vj^fc (11) 

d k , Vfc G /C (12) 



Remark: Suppose all the elements of the channel matrices 
are randomly and independently generated from continuous 
distributions and V^,V^,k G JC can be found to satisfy 
condition (11). Then, condition (12) will also be satisfied with 
probability 1. This is because the direct channel matrices Hl fefc l 
do not appear in condition (11). So the choice of transmit 
and receive filters V[ fc ],Tj[ fc ],fc G JC to satisfy (11) does not 
depend on the direct channel matrices H^- kk \ Since Hl fefe ] is 
independent of V^ k \XJ^ and all its elements are randomly 
generated from a continuous distribution (i.e. it lacks any 
special structure), the product matrix Tj[ fc l T H[ fcfc ] has full 
rank with probability 1. Thus, for random MIMO channels 
without time-extensions, if (11) can be satisfied then (12) is 
automatically satisfied almost surely as well. However, if time- 
extensions are considered then the channel matrices may have 
a block diagonal structure and (12) cannot be taken for granted. 

The solution to the feasibility problem is not known in 
general. In other words, given a set of randomly gener- 
ated channel matrices and a degree-of-freedom allocation 
(dW.dPl, • • • , d,W), it is not known if one can almost surely 
find transmit and receive filters that will satisfy the feasibility 
conditions. The distributed interference algorithm developed 
in this paper will be useful in (numerically) solving this open 
problem. 



B. Reciprocity of Alignment 

An interesting observation from the problem formulation 
above is the duality relationship between interference align- 
ment on a given interference channel and its reciprocal channel 
obtained fry_switching the direction of communication. Specif- 
ically, let V [ fc l , U denote the transmit precoding filters and 
the receive interference suppression filters on the reciprocal 
channel. The feasibility conditions on the reciprocal channel 



are: 



such that 



yW : N [k] x d [k] , U [felt U [fcl 

Ubltfjlife] v[ fe l 
rank(u[ fe ]tH[ fefe ]v[ fe ]) 



(13) 
(14) 



O.Vj^fc (15) 
d k , Vfc G K (16) 



Suppose we set = U^, tfM = V^. Then the feasibil- 
ity conditions on the reciprocal channel become identical to the 
original feasibility conditions. Thus, the following observation 
can be made: 

Reciprocity of Alignment: Since the feasibility condi- 
tions are identical, if the degrees of freedom allocation 
(dW , dffl , • • • , d\ K ^ ) is feasible on the original interference 
network then it is also feasible on the reciprocal network (and 
vice versa). Interference alignment on the reciprocal inter- 
ference network is simply achieved by choosing the transmit 
filters and the receive filters on the reciprocal channel as 
the receive filters and the transmit filters (respectively) of the 
original channel. 

Reciprocity of alignment is a key property used for dis- 
tributed interference alignment algorithms, described in the 
next section. 

V. Distributed algorithm for interference 

ALIGNMENT 

In this section we construct distributed interference align- 
ment algorithms for the interference channel with multiple- 
antenna nodes and no symbol extensions. Continuing with 
the system model of Section IV, this implies that the rele- 
vant interference alignment feasibility condition is (11) while 
(12) is automatically satisfied. Basically, (11) requires that 
at each receiver, all interference is suppressed, leaving as 
many interference-free dimensions as the degrees of freedom 
allocated to that receiver. 

Since we are interested in distributed algorithms, we start 
with arbitrary transmit and receive filters V[ fe ],U[ fc ] and itera- 
tively update these filters to approach interference alignment. 
The quality of alignment is measured by the power in the leak- 
age interference at each receiver, i.e. the interference power 
remaining in the received signal after the receive interference 
suppression filter is applied. The goal is to achieve interference 
alignment by progressively reducing the leakage interference. 
If interference alignment is feasible then eventually leakage 
interference will be zero. 

The total interference leakage at receiver k due to all 
undesired transmitters (j ^ k) is given by: 

I [k * ] = Tr ru [fclt Q [fel U [fcl j (17) 



where 



qW= ^H^V^V^H^t (18) 
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Reverse communication direction 



Fig. 2. Pictorial representation of the iterative interference alignment algorithm where the receive directions are optimized to minimize interference power 
at the receivers. Each link (arrow) represents a MIMO channel. The transmit power per node is P in both directions. 



is the interference covariance matrix at receiver k. 

Similarly, in the reciprocal network, the total interference 
leakage at receiver j due to all undesired transmitters (k ^ j) 
is given by: 



Tr 



Ub'ltQblubll (19) 



where 



K 



Q w = E 

fc=l,fe#j 



p[k] 
~dW 



Hbk] y[fc] y[fc]t jj[jfc]t (20) 



is the interference covariance matrix at receiver j. 

The iterative algorithm alternates between the original and 
reciprocal networks. Within each network only the receivers 
update their interference suppression filters to minimize their 
total leakage interference. 

Step I: In the original network, each receiver solves the 
following optimization problem. 



mm 

\j[k].]\[[k] X( [[k] Tj[ fc lUMt=I 



j[k* 



(21) 



In other words, receiver k chooses its interference suppression 
filter U^l to minimize the leakage interference due to all 
undesired transmitters. The d^ dimensional received signal 
subspace that contains the least interference is the space 
spanned by the eigenvectors corresponding to the smallest 
eigenvalues of the interference covariance matrix Q^. Thus, 
the d^ columns of are given by: 

Ufj = i/ d [QW], d=l, (22) 

where Vd[A\ is the eigenvector corresponding to the d th 
smallest eigenvalue of A. 

Step II: The second step is identical to the first step, but 
performed in the reciprocal network. Consider the reciprocal 



network obtained by reversing the roles of the transmitters and 
the receivers. The transmit precoding matrices in the recip- 
rocal network, V^, are the receive interference suppression 
matrices from the original network that were determined 
in Step I. Each receiver in the reciprocal network solves the 
following optimization problem. 



min_ _ I 



(23) 



Similar to Step I, the d^'l columns of U b'l are given by: 

U^] = ^[Q^], d=l, •■■>•] (24) 

The receive interference suppression filters in the reciprocal 
network are then used as the transmit precoding matrices 
in the original network, and the algorithm returns to Step 
I. The iterations continue in this manner until the algorithm 
converges. 

The iterative procedure is summarized in Algorithm 1. A 
pictorial representation is shown in Fig. 2. 

A. Proof of Convergence 

We now show that the algorithm must converge. The proof 
also highlights the intuition behind the algorithm. 

We define a metric called the weighted leakage interference 
(WLI) as: 

K <p[fc] 



III) 



K 

fc=l j=l,jjtk 

K K 



E 



.jim 



- E E 



pw pu] 



■Tr 



u [ fe ]tH[ fej ]V [j ' 1 V [j ' lt H to ' lt U [,cl 



k=lj=l,j^k 

We show that each step in the algorithm reduces the value 
of WLI. Since WLI is bounded below by zero, this implies 
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Algorithm 1 Iterative interference alignment 
1: Start with arbitrary precoding matrices V^'l : x 




Therefore the value of U^l computed in Step 7 to minimize 
I [■ J *l also minimizes Iw^. Since /„, = J2k=i ^\ it is easily 
seen that Step 7 can also only reduce the value of I w . Since 
the value of /„, is monotonically reduced after every iteration, 
convergence of the algorithm is guaranteed. 

Remark: While the algorithm minimizes leakage interfer- 
ence at every iteration and is guaranteed to converge, con- 
vergence to global minimum is not guaranteed due to the 
non-convex nature of the interference optimization problem. 
Numerical results for the performance of the algorithm are 
presented in the next section. 

The following observations summarize the intuition behind 
the iterative algorithm. 

1) Dimensions along which a receiver sees the least inter- 
ference from other nodes are also the dimensions along 
which it causes the least interference to other nodes in 
the reciprocal network where it functions as a transmitter. 

2) The weighted leakage interference is unchanged in the 
original and reciprocal networks if the transmit and 
receive filters are switched. 



2: Begin iteration 

3: Compute interference covariance matrix at the receivers: 
Q [k] = ^H^'lV^V^tH^'lt 

4: Compute the interference suppression matrix at each re- 
ceiver: 

U™=MQ l % d =1, •••,#] 

5: Reverse the communication direction and set = U^. 
6: Compute interference covariance matrix at the new re- 
ceivers: 

1: Compute the interference suppression matrix at each re- 
ceiver: 

8: Reverse the communication direction and set Vl fc l = U^L 
9: Continue till convergence. 



that the algorithm must converge. Note that an interference 
alignment solution corresponds to WLI= 0. 
The WLI associated with receiver fc is 

*E[k] K phi r 



d [k] ^ d [j] 

p[k] 



■Tr 



p[k] 

~dW 



j[k*] 



Therefore the value of computed in Step 4 to minimize 
a i so 

minimizes Iw K Since I w — X^fe=i ^ ^> we have 



min I w = 

UW.UPl,- ,UM 



K 

min &* 



fe=l 



K 

E 

fc=i 



min I [k * 
uw 



1^ d [k] 

k—1 



min 7^ 

ui'i 



In other words, given the values of , j £ {1, 2, • • • , K }, 
Step 4 minimizes the value of I w over all possible choices of 
UM , k e {1, 2, • • • , K}. In particular, Step 4 can only reduce 
the value of I w . 

The weighted leakage interference associated with transmit- 
ter j is 



P [j] * p[k] 



dbl f-f dM 

fe=l 



p[i] 



if 

E 



p[k] 



M dM 
fe=i 



■Tr 



V [fc] t h \i fe ] t u b'l U yi f H [ J ' fe ] V [fe] 



Max-SINR Algorithm 

The algorithm presented above seeks perfect interference 
alignment. In particular it seeks to create an interference- 
free subspace of the required number of dimensions, that 
is designated as the desired signal subspace. However, note 
that interference alignment makes no attempt to maximize 
the desired signal power within the desired signal subspace. 
In fact the algorithm described above does not depend at 
all on the direct channels H[ fefe l through which the desired 
signal arrives at the intended receiver. Therefore, while the 
interference is eliminated within the desired space, no coherent 
combining gain (array gain) for the desired signal is obtained 
with interference alignment. While this is optimal as all 
signal powers approach infinity, it is not optimal in general 

5t intermediate SNR values. Therefore other algorithms may 
e designed which will perform better than the interference 
lignment algorithm at intermediate SNR values. 
In this section we consider one such natural extension of 
the interference alignment algorithm where the receive filters 
\j[k] and Tj[k] 

are chosen to maximize SINR at the receivers 
instead of only minimizing the leakage interference. While 
there is no loss of generality in assuming orthogonal precoding 
vectors for the streams sent from the same transmitter as far 
as interference alignment is concerned, orthogonal precoding 
vectors are in general suboptimal for SINR optimization. We 
therefore no longer assume that the columns of (the 
transmit precoding vectors) are mutually orthogonal. We also 
identify the columns of to be the specific combining 
vectors for the corresponding desired data stream, so that 
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they are not necessarily orthogonal either. With these modified 
definitions, the SINR of the I th stream of the k th receiver is 



TJ [k] t TJ [fefe] V [fe] V [fe] t TJ [fefe] t Tj[fc] p[k] 

" -kl *l 



u^ ]t B[*qug ] 



SINR/,, 
where 

B[H1 = ESrE Hto ' Iv W Hto1t 



(25) 



k p[j] &l 



-^H^V^vifH^t+i^ (26) 

\k] 

The unit vector LP ; that maximizes SINRfe; is given by 

1 H -" " V - 

(27) 



B 



ki] 



H [fcfc]y[fe] 



(B[ fei ]) _1 H[ fcfe ]V^ ||' 
The steps of the iteration are given in Algorithm 2. 

Algorithm 2 Max-SINR algorithm 

1: Start with any V[ fc l : x <#1, columns of V^l are 

linearly independent unit vectors. 
2: Begin iteration 

3: Compute interference plus noise covariance matrix for 
g[fci] f or s tream I at receiver fc according to (26), V fc G 
{1,2,-.. ,x},;g {1,2,--- ,#]}■ 

4: Calculate receive combining vectors at receiver fc ac- 
cording to (27), VfcG {1,2,- ■• ,^},/G {1,2,- ■• ,#]}. 

5: Reverse the communication direction and use the receive 
combining vectors as precoding vectors: = U™, 

6: In the reciprocal network, compute interference plus-noise 
covariance matrix B ^ for stream I at receiver k, V k G 
{1,2 1 - ) 4!6{1,2,...,#]}. 

ffel 

7: Calculate receive combining vectors U 1 , , V k £ 

{1,2,... ,X},^G {1,2,-.. ,#]}. 

8: Reverse the communication direction and use the receive 
combining vectors as precoding vectors: V^'l = U^, 
Vfce {1,2,-.. ,K}. 

9: Repeat until convergence. 



VI. Performance Results and Applications 

Consider the 3 user interference channel where each node 
is equipped with 2 antennas and all channel coefficients are 
i.i.d. zero mean unit variance circularly symmetric complex 
Gaussian. As shown in Fig. 1 with interference alignment each 
user achieves 1 degree of freedom. In Fig. 3 we compare the 
performance of the iterative interference alignment algorithm 
for this channel with the theoretical interference alignment 
which assumes global channel knowledge. It can be seen 
that the algorithm performs very close to the theoretical case 
and more importantly it provides significant benefits over the 
orthogonal case. The sum rate for the orthogonal scheme is 
calculated assuming equal time sharing for the users, and 
with power 3P per node. The modified interference alignment 
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Fig. 3. Performance of the decentralized interference alignment algorithm 
for the three user two antenna case. 



scheme where the receive-combining vectors maximize the 
SINR provides a considerable performance gain at low and 
intermediate P. As expected, it achieves the same performance 
of interference alignment at high P. We also plot the perfor- 
mance of interference avoidance algorithm which is a selfish 
approach. The isotropic transmission case refers to the case 
where each transmitter sends 2 streams of equal power without 
regard to the channel information. In the following, we explore 
some of the applications of the iterative interference alignment 
algorithm. 

A. Feasibility of Interference Alignment 

While the iterative algorithm is useful for circumventing 
the need for global channel knowledge, it can also be used 
to check theoretical feasibility of interference alignment for 
a given number of streams per user. Let (g?W , , • • ■ d^ ) 
denote the number of transmit streams of the users. For perfect 
interference alignment ^2j=i [Q^] = at receiver k where 
Aj[A] denotes the jth smallest eigenvalue of A. Note that 

J2j=i ^jlQ^] indicates the interference power in the desired 
signal space. 

Using the algorithm, we plot in Fig. 4, the percentage of 
interference in the desired signal space versus the total number 
of transmit streams in the network. The fraction of interference 
in the desired signal space of receiver k is defined as 



Pk 



[Q 



Ml 



(28) 



Tr[QW] ' 

When interference alignment is feasible the fraction of inter- 
ference in desired signal space will be zero (within numerical 
errors). Fig. 4 suggests that interference alignment is feasible 
on the four user interference channel with 5 antennas at 
each node when each transmitter sends two streams. With 
increase in the number of streams the interference in desired 
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Fig. 4. Percentage of interference power in desired signal space as a function 
of the total number of data streams in the network. 



signal space increases which is an indication that interference 
alignment is not possible. Although the upperbound on the 
degrees of freedom for this network is 10, approaching the 
upperbound may entail channel extensions. The plot suggests 
that 8 degrees of freedom (d^ = dfi = d^ = efM = 2) 
can be achieved without channel extension. Similarly for the 
4 antenna case, the plot indicates that interference alignment is 
possible for up to a total of 6 streams in the 4 user interference 
network with only 4 antennas at each node. 

B. Networks with single antenna nodes 

MIMO nodes are not necessary in order to achieve inter- 
ference alignment in wireless networks. Interference can be 
aligned even in networks with single antenna nodes through 
channel extension in frequency or time as long as the channel 
is varying across frequency or time [1]. However, one caveat 
with this approach is the need for long symbol extensions. 
Aligning interference with a large number of beams over a 
large number of dimensions can be especially challenging 
for iterative algorithms due to large dimensionality of the 
optimization space and the inherently non-convex nature of 
the problem. Moreover, symbol extensions over orthogonal 
dimensions produce structured (diagonal or block diagonal) 
matrices for which both conditions (12) and (11) are non- 
trivial. Due to these difficulties, it is preferred if interference 
alignment can be accomplished without symbol extensions or 
with limited symbol extensions. 

In this section, we give an example of how long symbol 
extensions may be avoided by the use of relays. Note that [16] 
has shown that relays cannot increase the degrees of freedom 
for time-varying wireless networks. However, as we show in 
this section, relays can be very useful by reducing the size of 
the signalling space over which interference alignment can be 
accomplished. The key idea is to employ relays to create a 
virtual MIMO system. Consider an interference relay channel 
with three sources, three destinations and a half-duplex relay 
(node 0) as shown in Fig. 5. Recall that in the absence of 



relays this network is shown to approach the upperbound 
of 3/2 degrees of freedom per orthogonal dimension in the 
asymptotic limit of infinitely long symbol extensions [1]. 
However, we show that with relays only a two time slots 
are required to achieve the outerbound, i.e. 3/2 degrees of 
freedom. 

Consider the following two slot protocol. In the first slot, 
the relay is silent and the received signal at destination j is 
given by 

3 

ybl(i) =J2h m (l)x®(l) j = 1,2,3 

»=i 

The received signal at the relay can be expressed as 

3 

yM(l) = Y / hm ( 1 ) xll] ( 1 ) + zla] ( 1 ) 

2 = 1 

In the second slot, the relay transmits a scaled version of its 
received symbol while source i transmits x^(2). The received 
signal at the j th destination node in the second slot (j = 

1,2,3) 

3 

y [l] (2) = hM(2)X®(2) + h^(2)/3y^(l) + z^(2) 
j=i 

Let = [y^(l) y^{2)] T and X® = [x®(l) x^(2)] T . 
In vector form, the received signal at destination j can be 
expressed as 

3 



(29) 



where 



and 
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f3h lj0] {2)h [0l] (l) /i M (2) 
zbl(i) 

f3W (2)zW(l) + z^(2) _ 



Thus, over two time slots, the relay network reduces to a 
three user MIMO interference channel with (non-diagonal) 
structure on the channel matrix. Since the channel matrix 
is non-diagonal (unlike symbol extensions in the absence 
of a relay) it is easy to verify that when the channels are 
random and independent of each other, a multiplexing gain 
of | is achieved with probability 1. The advantage of this 
scheme is that it requires only two time-slots to achieve 3/2 
degrees of freedom per time-slot, whereas in the absence of the 
relay, infinitely many time-slots are used in [1]. Since fewer 
dimensions are needed, iterative algorithms work better (faster 
convergence) in this setting. 

In Fig. 6 we plot the performance of the interference 
alignment schemes for the case of time extensions. For the 
interference alignment scheme with relay, the transmit power 
per node is P. That is, the three transmitters and the relay have 
a total power of 4P. The transmit power for the orthogonal 
scheme where only transmitter is active at a time is 4P. 
The interference alignment scheme without relay achieves a 
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VII. Conclusion 

Distributed interference alignment algorithms are investi- 
gated. Interference alignment is found to be achievable through 
iterative algorithms based on network reciprocity and the 
"minimize interference to others" approach. Numerical com- 
parisons to orthogonal schemes, simultaneous transmission 
schemes and selfish interference avoidance schemes show that 
the benefits of distributed interference alignment algorithm 
are significant and close to the theoretical predictions. As 
mitigating interference is the fundamental problem of wireless 
networks, the 'do no harm' approach based algorithms have 
enormous applications in wireless networks. 
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Fig. 6. Performance of the interference alignment schemes with time slot 
extension for the three user interference channel. 



multiplexing gain of 4 in 3 time slots [1]. Here transmitter 1 
sends two streams while transmitters 2 and 3 send one stream 
each. The transmit power per stream is P adding up to 4P for 
the scheme. We can also increase the number of time slots for 
channel extension to improve the multiplexing gain. However 
it requires very high P to outperform the orthogonal scheme. 
It can be seen that adding an extra relay helps in achieving the 
multiplexing gain of | in two slots. Further the performance 
improves at low P as well. 

It must be stressed that the main idea behind the interference 
alignment scheme with the relay is to show that there are ben- 
efits in employing relays when the nodes do not have multiple 
antennas, especially for iterative algorithms. The scheme can 
be further improved by optimally allocating power to the relay 
and employing Algorithm 2 to improve its performance at low 
P. 
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